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We study the statistical properties of one-dimensional directed polymers in a short-range random 
potential by mapping the replicated problem to a many body quantum boson system with attractive 
interactions. We find the full set of eigenvalues and eigenfunctions of the many-body system and 
perform the summation over the entire spectrum of excited states. The analytic continuation of 
the obtained exact expression for the replica partition function from integer to non-integer replica 
parameter A'^ turns out to be ambiguous. Performing the analytic continuation simply by assuming 
that the parameter A'^ can take arbitrary complex values, and going to the thermodynamic limit of 
the original directed polymer problem, we obtain the explicit universal expression for the probability 
distribution function of free energy fluctuations. 
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I. INTRODUCTION 



Directed polymers in a quenched random potential have been the subject of intense investigations during the past 
two decades [i| . Diverse physical systems such as domain walls in magnetic films [1] , vortices in superconductors 0j , 
wetting fronts on planar systems or Burgers turbulence Q can be mapped to this model, which exhibits numerous 
non-trivial features deriving from the interplay between elasticity and disorder. The best understanding has been 
achieved for a string confined to a plane. In this case we deal with an elastic string directed along the r-axis within 
an interval [0, L]. Randomness enters the problem through a disorder potential y[0(r),T], which competes against 
the elastic energy. The problem is defined by the Hamiltonian 



H[cI>{t),V] 



dT[\[dr 



v\ 



(1) 



where in the simplest case the disorder potential r] is Gaussian distributed with a zero mean r) — and 
the (5-correlations: 



V{(j), T)V{(j)', t') = u5{t - r')5(0 - 0') 



(2) 
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Here the parameter u describes the strength of the disorder. Historically, the problem of central interest was the 
scaling behavior of the polymer mean squared displacement which in the thermodynamic limit {L — >■ oo) is believed 
to have a universal scaling form 



(02)(L)cxL2C (3) 



(where (...) and (...) denote the thermal and the disorder averages), with C, the so-called wandering exponent. More 
general problem for all directed polymer systems of the type, Eq.([T|), is the statistical properties of their free energy 
fluctuations. Besides the usual extensive (linear in L) self-averaging part f^L (where /o is the linear free energy 
density), the total free energy F of such systems contains disorder dependent fluctuating contribution which is 
characterized by non-trivial scaling in L. It is generally believed (at list in the systems with short-ranged correlations) 
that in the limit of large L the typical value of the free energy fluctuations scales with L as 

F cx L"', (4) 

i.e. they are characterized by a single universal exponent uj < 1. In other words, in the limit of large L the total 
(random) free energy of the system can be represented as 

F = /oL + (5) 

where / is the random quantity which in the thermodynamic limit L ^ oo is described by a non-trivial universal 
distribution function P*{f)- The derivation of this function for the system with (5-correlated random potential, 
Eqs.([T])-(l2|) is the central issue of the present work. 

One can easily note that the above two exponents C a-nd w are not independent. Indeed, since the free energy 
fluctuation F ^ L'^ can be estimated by the typical value of the elastic energy, F cx (j)^ /L, where, according to Eq.Q, 
the typical deviation (f> ^ , one finds a; = 2^ — 1. 

It is generally believed that for all short-range correlated disorder potentials, the free energy fluctuations exponent 
assumes a universal value uj — 1/3. Numerical studies [6| as well as the solution via mapping to the Burgers equation 
Q confirm this conjecture. One arrives to the same conclusion studyingscaling properties of the free energy by 
mapping the replicated problem to N particle quantum bosons system (8| and using the Bethe-Anzats solution. 
However, in this latter case, the resulting distribution function P*{f) exhibits severe pathologies such as the vanishing 
of its second moment, which assumes that the distribution function is not positively defined. 

Let us consider this point in more detail. For the string with the zero boundary conditions at r = and at r = L 
the partition function of a given sample is 

z[v] = / v[(j){T)] (6) 

J0(O)=O 

where /3 = 1/T denotes the inverse temperature. On the other hand, the partition function is related to the total free 
energy F[V] via 

Z[V]^eM~l3F[V]) (7) 

The free energy F[V] is defined for a specific realization of the random potential V and thus represent a random 
variable. Let us take the A^-th power of both sides of Eq.(I7]) and perform the averaging over the random potential V: 



Z[N,L]= [eM-mF[V])] (8) 
The quantity in the l.h.s of the above equation 



Z[N,L] = Z^[V] (9) 

is called the replica partition function, and it is defined originally for an arbitrary integer parameter N. Substituting 
F ^ foL + fL'^, into Eq.® and redefining 

Z[N, L] ^ Z[N, L] Q-P'^^o^ (10) 

we get 



^[iV,L]= exp(-/?iVL"/) (11) 
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The averaging in the r.h.s of the above equation can be represented in terms of the distribution function -Pl(/) (which 
depends on the system size L). In this way we arrive to the following general relation between the replica partition 
function Z[N,L] and the distribution function of the free energy fluctuations PL{f)- 

/ + 00 
dfPL{f)e-^^'^^f (12) 
-oo 

The above equation is the bilateral Laplace transform of the function PL{f), and at least formally it allows to restore 
this function in terms of the replica partition function Z[N,L]. In order to do so we have to compute Z[N,L\ for 
an arbitrary integer N and then perform analytical continuation of this function from integer to arbitrary complex 
values of N . Introducing a new complex variable 

S = PNL"^ (13) 

and denoting 

Z[^,L]^Zl{^ (14) 
we could reconstruct the distribution function PL{f) via the inverse Laplace transform 

PlU) = / —ZL{s)e'f, (15) 

J_ioo 27^^ 

Finally, provided there are exist a finite thermodynamic limit function 

lim Zl{S) = Z,{S) (16) 

we can find the distribution function 

r+ioo i~ 

which would describe the statistics of the rescaled free energy fluctuations / in the infinite system. The above 
equation defining P*{f) contains no parameters and hence is expected to be universal. Therefore according to the 
relation s = /3NL'^ we see that in the thermodynamic limit the relevant values of the original replica parameter are 

N - L-'^ (18) 

This explains why the two limits L oo and iV — > do not commute Q, and the approximation of the replica 
partition function through the ground state wave function fails (see also Ref. [13] )• In Kardar's original solution Q, 
after mapping the replicated problem to interacting quantum bosons, one arrives at the replica partition function for 
positive integer parameters iV > 1. Assuming a large L — > oo limit, one is tempted to approximate the result by 
the ground state contribution only, as for any N > 1 the contributions of excited states are exponentially small for 
i — ^ oo. However, in the analytic continuation for arbitrary complex N the contributions which are exponentially 
small at positive integer > 1 can become essential in the region — 0, which according to Eq. (IT7|) defines the 
function P*(/) Thus, it is the neglection of the excited states which is the origin of non-physical nature of the obtained 
solution. In other words, for the proper analytic continuation of the replica solution to the region A^ — ^ 0, first one 
has to calculate the replica partition function Z[N,L] exactly for arbitrary integer A^, and only after that one can 
take the thermodynamic limit L — > oo while keeping the value of the parameter s = (3NL'^ finite. 

In the present paper, we report the results of the calculation of the replica partition function Z(N, L) for arbitrary 
integer A^ which in terms of the Bethe-Ansatz solution for quantum bosons with attractive (5-interactions (Sections 
II - III) involves the summation over the entire spectrum of exited states (Section IV). Unfortunately the analytic 
continuation of the obtained exact expression Z(N,L) from integer to non-integer A turns out to be ambiguous, 
since our replica partition function growth as exp(A^'^) at large A' (similar problem of the analytic continuation to 
the region N ^ one faces in the replica theory of the mean-field spin-glasses where the replica partition function 
growth as exp(A'^) [lH). Performing a kind of a "replica symmetric" analytic continuation, i.e. just assuming that 
originally integer-value parameter A^ can take arbitrary complex values and taking the thermodynamic limit L — >■ oo 
in Z{N, L) (Section V) allows us to compute an inverse Laplace transformation, cf. Eq. (I17|) . which provides us with 
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the explicit expression for the distribution function of the free-energy fluctuations (Section VI, Eq. (|102l) l. Although 
up to the present moment, we have not uncovered any unphysical properties in the obtained probability function 
P*{f), this solution could be considered as distant analog of the "replica symmetric approximation" in the mean-field 
spin-glasses. In particular, it should be noted that our result is different from the Tracy- Widom distribution fl^ . 
which describes the statistics of fluctuations in various statistical systems jl3r,18] which are widely believed to belong 
to the same universality class as the present model (l9l - l2l| (for further discussion of this issue see Section VII) . 

Various technical aspects of the calculations are moved to the Appendices. In particular, in Appendices A, B, and 
C, we analyze the structure and properties of A^-particle wave functions of one-dimensional quantum bosons, both 
with repulsive and with attractive interactions. 



II. MAPPING TO QUANTUM BOSONS 

Explicitly, the replica partition function, Eq.Q, of the system described by the Hamiltonian, Eq.([T]), is 



N 



Z{N,L) = l[ 



0a(i)=O 



'D(j)a{T) exp 



0a(O)=O 



N 



(19) 



Since it is assumed that the random potential V\4>tt\ has the Gaussian distribution the disorder average (...) in the 
above equation is very simple: 



exp 



-P / dT^l/[0„(r),T] = exp ^ / / drdr' V V [0jr), T]F[0b(T'), r'] 



(20) 



Using Eq.© we have: 



z(iv,L)=n 



Vc/fair) exp 



-/3 / dT[J2[drMr)]' - Pu J2 S[Mr)-Mr)]} 



a,b=l 



(21) 



It should be noted that the second term in the exponential of the above equation contain formally divergent contribu- 
tions proportional to (5(0) (due to the terms with a = b). In fact, this is just an indication that the continuous model, 
Eqs.((T|)-(l2]) is ill defined as short distances and requires proper lattice regularization. Of course, the corresponding 
lattice model would contain no divergences, and the terms with a = 6 in the exponential of the corresponding replica 
partition function would produce irrelevant constant ^Lf3'^uNS{0) (where the lattice version of (5(0) has a finite value). 
Since the lattice regularization has no impact on the continuous long distance properties of the considered system this 
term will be just omitted in our further study. 

Introducing the A''-component scalar field replica Hamiltonian 



Hn[4>] = - / driY.[9rUr)]^ - l3uY.S[Mr) - Mr)] 

\a=l a^b 



(22) 



for the replica partition function, Ea.(|2ip. we obtain the standard expression 



ZiN,L) = l[ 



(23) 



where </> = {(/)i,...,(/)jv}. According to the above definition this partition function describe the statistics of N 
(5-interacting (attracting) trajectories 4)a{T) all starting (at t — 0) and ending (at t — L) dX zero: (t>a{^) — 4>a{L) = 
In order to map the problem to one-dimensional quantum bosons, let us introduce more general object 



*(x;i) = n 



0a(O)=O 



v^t) e 



(24) 



which describes A'' trajectories (pair) all starting at zero ((^a(0) = 0), but ending at r = i in arbitrary given points 
{xi, a; at}. One can easily show that instead of using the path integral, ^'(x; t) may be obtained as the solution of 
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the linear differential equation 

N 1 ^ 

9t*(x;0 = J2dl^i^;t) + ^P^uJ2Sixa-Xb)^iK;t) (25) 

^ a=l a^b 

with the initial condition 

^{^■,0)=U^^,6{xa) (26) 
One can easily see that Ea. (|25l) is the imaginary-time Schrodinger equation 

-5t*(x;t) =iJ^'(x;i) (27) 

with the Hamiltonian 

^ = E ^1 - ^l3'uJ2S{xa - xb) (28) 

a=l 



which describes N bose-particles of mass /3 interacting via the attractive two-body potential —l3'^uS{x). The original 
replica partition function, Eg. (1231) . then is obtained via a particular choice of the final-point coordinates, 

Z(iV,i) = vI/(0;L). (29) 

The standard general strategy of the further calculations is in the following. Let us denote the eigenfunctions of the 
Hamiltonian Eq. (f28| by ^'§(x) where the index ^ (which can be both integer and continuous) labels the eigenstates. 
Provided the wave functions ^'^(x) constitute the orthonormal and complete set, the time dependent solution of the 
equation ((25| with the initial conditions, Ea. (l26l) . is given by 



*(x;t) = J2 *€W e--^(«)* (30) 

where E{^) denotes the energy of the ^-th eigenstate: 

ffvl,^(x) = i?(e)*c(x) (31) 

Then, according to Eg. d^Hl) . the replica partition function Z{N,L) of the original polymer system is obtained just by 
the summation over all eigenstates of the quantum Hamiltonian (1281) : 



Z{N,L) = J2 |*c(0)|^ e-^(«)^ (32) 

Thus, the crucial point of the present approach is finding the eigenfunctions and the energy spectrum of the Hamil- 
tonian (pS)) . which is the topic of the next section. 

III. EIGENSTATES OF THE ONE-DIMENSIONAL QUANTUM BOSONS SYSTEM 

A. Repulsive bosons 

The eigenfunctions of one-dimensional ^-interacting repulsive (u < 0) quantum bosons, Eq. (j28|) . have been derived by 



Lieb and Linigcr in 1963 |22| . An eigenstate of this system is characterized by iV continuous momenta {qi, qn} = q 



with the wave function (sec Appendix A) 

. N s N 

vl/^^)(x) = C(^)(q) ^ (-1)1^1 n - 9pJ + ^^ssn{xa - x,)] J exp[z ^ (33) 
where we have introduced the notation 

K = (34) 
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The summation in Ea. (p3)) goes over all permutations P of the A'' momenta {qi, ...,(7jv} over N particles located at 
{xi, ...,X7v}, and [P] denotes the parity of the permutation. The normalization constant C^^^(q) is 



C(^)(q) = , (35) 



and the associated energy £'Ar(q) is 



1 ^ 

EN{ci)^—Y^ql (36) 



2/3 

^ a=l 

A useful alternative representations of these wave functions are 

^ N s AT 

*W(x) ^ CW(q)^(-l)[^lf J][-z(5,„-a,J+*A^sgn(x,-x,)]j exp[z^gp„x,] 



(37) 



P ^a<b 
N 



= &''\ci)(\{[~i{d^^-d^,,)+iKSgn{xa-Xb)\ \ det[exp(iqx)] (38) 

where the symbol exp(iqx) denotes the N x N matrix with the elements exp{iqaXb) {a,b — 1,...,N). Here, by 
definition, the differential operators dx^ act only on the exponential terms and noton the signum functions sgn(a::a — Xb). 

One can easily see that the above wave functions ^'q^'' (x) represent a set of plane waves in any sector of the type 
Xai < Xa2 < ■■■ < Xaj^ with a finite jump (equal to k) of the derivatives at all "boundary" points Xai — Xa^. These 
functions are symmetric with respect to any permutation of the particle coordinate {xi, ...,xn} and antisymmetric 
with respect to permutations of the momenta {qi, qjsr}. It can be proven that the wave functions, Eq. p3p . are 
orthonormal and constitute the complete set (see Appendix A; a detailed discussion one can found e.g. in Refs. 

[l^ll^l)- Specifically, for any two functions 4'q^''(x) and ^'^^■'(x) considered in the sectors qi < q2 < ■■■ < qN and 
q'l < q'2 < ... < q'jsj orthonormality implies that 

/+00 
dx^...dxN <^)*(x)vI/W(x) = (2^)^<5(g; - q^)5{q'2 - q2)...5[q'^ - g^) (39) 
-00 

Similarly, for any two functions ^q^^(x) and ^q^^(x') considered in the sectors xi < X2 < •■• < xn and x'-^ < x'2 < 
... < x'^, completeness implies that 

dq^...dqM *W-(x)vI/W(x') = (2^)^5(xi - x'^)5{x2 - x'2)...S{xN - x'j,) (40) 

-oo 

Thus, for a repulsive interaction {u < 0) the time dependent solution of the differential equation ((25)) (with the starting 
condition, Eq. ipS)) ) would be sufficiently simple: 

^^"'(x;0 = £^ ^^|^*r*W*r(0)exp[-±f;,^] (41) 

Unfortunately, from the point of view of the replica theory of disordered polymers the repulsive bosons make no 
physical meaning since the parameter u (according to its definition, Eq.([5])) is positively defined. 



B. Attractive bosons 



The situation with attractive {u > 0) bosons is more complicated. One can easily prove that (irrespective of the 

sign of the parameter u) the functions \I'q^''(x), Ea. (l33|) . are othonormal eigenfunctions of the Hamiltonian, Eg.pS]). 
However, unlike the repulsive case, this set of functions is not complete. In other words, for u > the completeness 
conditions, Eqs. lHO)) . are not satisfied. Physically this indicates that besides the continuous spectrum (or free particles) 
states, our system must have another types of eigenstates, in which the particles are bound into localized clusters. 
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The spectrum and some properties of the eigenfunctions for attractive one-dimensional quantum bosons have been 
derived by McGuire [1^ and by Yang [2^ (see also Ref. [13, Ull). However, since attractive bosons do not have a 
proper thermodynamic limit (in the number of particles N oo) due to the scaling oc —N^, the interest in this 
system has been rather limited. 

We first consider the ground state wave function ^^^■'(x) in which all N particles are bound into one cluster with 
the free center of mass motion controlled by the momentum q (see Appendix B): 



*W(x) = C^^\q) exp 



N 

iq'^Xa 



1 



N 



(42) 



where 



kN 



(43) 



is the normalization constant defined by the orthonormahty condition 



+ 00 



*W(x)*l'^*(x) = (27r)(5(g - g') 



(44) 



The energy of this state is 



Eiiq;N) 



TV 

2^^ 24/3 



{N^ - N) 



(45) 



On the other hand, one can also easily prove that the above ground state wave function, Ea. (|42)) . can be represented 
in the form similar to the free particle structure, Ea (j33|) . by introducing (discrete) imaginary parts for the momenta 
qa- Indeed, due to the symmetry of this function with respect to permutations of {a;i,a;2, ...,a;jv} it is sufficient to 
consider it in the sector a;i<X2<...<a;jv. Defining 



qa = q--K{N+l-2a) 



(46) 



and substituting these momenta into the general expression for the wave function, Ea. p3p one easily recovers Eq. (j42p 
(see Appendix B for details) Also, substituting Ea. (l46t into the general expression for the energy spectrum. Eg . ((36)) 
one can also recover, Ea. (H5t . 

Now, using the above scheme, one can construct the eigenfunctions of a generic excited state. It consists of 
M (1 < M < N) "clusters" {ila} of bound particles, where a = 1,...,M labels a given cluster. Each cluster is 
characterized by the momentum qa of its center of mass motion, and by the number Ua of particles contained in it 



M 



(such that 
"vector" momenta 



N). Instead of N independent real momenta qa 



l,...,N) one introduces M complex 



qa 



qa 



K [ria + 1 — 2r) 



(47) 



where 



1, 2, Ua and 



M 

E 

a=l 



N 



(48) 



The corresponding wave function 'i'^n (xi, ...,xn) is characterized by M continuous parameters q — (qi, ...,qM) 
(which are the momenta of the center of mass motion of the clusters) and M integer parameters n = (rii, ...jUm) 
(which are the numbers of particles of each cluster) . Explicitly this wave function is given by Eq. ([55]) where N x N 

matrix exp(iqx) is now composed of N columns 



{ql, ql, 



qln ; 



ql 



ql 



ql2' 



q¥. 



q^' 



} 



(49) 



and N rows Xa- One can easily see that as in the case of repulsive bosons, this wave function is symmetric with respect 
to permutation of particles coordinates. However, for practical applications the general representation, Eq. (j38|) . is 
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not very convenient. Writing the determinant of the matrix exp(zqx) explicitly, after a few efforts in simple algebra 
one can derive more transparent structure of the wave function (see Appendix C for details). Assuming first that the 
position of particles are ordered, xi < X2 < ... < xjv, let us consider a permutation P of N momenta q", Ea. (Hn|) . 
over N particles Xa, so that a particle number a is attributed a momentum component q"l^^y The particles getting 
the momenta with the same a (having the same real part Qa) will be called belonging to a cluster 17^. For a given 
permutation P the particles belonging to the same cluster are numbered by the "internal" index r = 1, riQ,. From 
now on we have to take into account not all the permutations, but only those for which the internal particle number 
r{a) is the growing function of the particle number a in every cluster. Namely, let a cluster fla consists of the particles 
Xai , , ■ ■ ■ 7 Xa„^ (whcrc fli < 02 < • • • < a„^ and the positions of particles are ordered: Xa-^ < < ■ ■ ■ < Xa„^ ). 
Then, among all Ual "internal" permutations of the momenta components qf, , • ■ ■ , non-zero contribution is 
given only by the one in which r(aj+i) — r{ai) + 1. In this case the explicit form of the wave function (for M > 2) 
reeds (see Appendix C) 



*rj(x) = (-1)1^1 n 



a<b 



1a{a) - ■^K("a(a) + 1 - 2r(a)) ) - ( q^ib) - ^K(n-a(fc) + 1 - 2r(6)) ) - in 



X exp 



N 



N 



(50) 



where the product goes only over pairs of particles belonging to different clusters and the symbol X]p' nieans that the 
summation goes only over the permutations P in which the "internal" indices r(a) are ordered inside each cluster. 

For generic positions of the particles (beyond the sector xi < X2 < ■■■ < xn) the expression for the wave function 
reduces to 



JV 



a<b 
a{a)y^a{b) 



qa{a) - 9a(6) + X! Sgn(Xa " a^c) " ^ ^^^'^i^b ^ ^c) + IK Sgn{Xa - Xb) 



(i>) 



X exp 



M M 



a— 1 a^Oa 



a— 1 a^a'i^VLa 



(51) 



Here the summation goes only over the permutations in which "internal" indices r(a) in the clusters are ordered 
according to the spatial ordering of the particles belonging to these clusters. For example, let a cluster $7^ be composed 



of the particles {xa^ , 



} and the spatial positions of these particles are such that Xa-^ < Xa2 < ■ ■ ■ < Xa„ 



while the particles numbering 01,02, .. . , o„^ is now arbitrary. Then the ordering of the internal index r{a) of the 
permutations involved in Ea. (j5ip is such that r(ai+i) — r{ai) + 1 (i.e. it goes from the smallest Xa in the cluster to 
the largest one). 

One can easily see that the expression for "i'q^{:x.), Eg. (|5T|) . can also be rewritten in more compact form: 



*rJ(x) = C(^)^(-l 
p 



N 

n 

a<b 
a{a)^a{b) 



-i{dx^-dx^) +iKSgn{xa-Xb)] 



exp 



■ M 

*E«" E 

■ a=l aGHo 



M 
K \ ^ 



E 

Q — 1 a,a'Gfla 



(52) 



In order to describe the orthogonality of these wave functions, we have to specify their symmetry structure. First 
of all, like in the case of repulsion. Eg. (155)) . the wave function of unbound particles (the case M = N and ni = n2 = 
... = — 1) is fully antisymmetric with respect to any momenta permutation. For the generic case, 2 < M < N the 
situation is slightly more tricky as the symmetry with respect to the momenta permutations depends on the values 
of the corresponding integer parameters Ua- According to Eq.(|38|). the permutation of any two momenta qa^ and 
qa2 belonging to the clusters which have the same numbers of particles, produces the factor (— 1)" 

(this operation corresponds to the permutation of n columns of the matrix exp{iqaXb)) ■ Hence two eigenstates which 
differ one from another only by such momenta permutations could be called equivalent. On the other hand, the 
permutation of the momenta between two clusters with different numbers of particles reveals no specific symmetry 
at all. In other words, the wave functions with permuted momenta belonging to two clusters with different numbers 
of particles are just two different wave functions describing two different eigenstates. Thus for comparing the wave 
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functions described by the parameters {qa,na) {a = 1, M) (in particular for the study of their orthogonahty) it 
is crucial to specify "subsets" of equal n's. Namely, a generic eigenstate (q, n) with M clusters could be specified in 
terms of the following set of parameters: 

(q,n) = {iqi,mi),...,{qs,,mi);{qs,+i,m2),...,{qsi+s2:m2); .... ; ((7si+...+st_i+i, Wfe), , Wfe)} (53) 

^ V ' V ' . ' 

Si S2 Sk 

where Si {i — 1, ...fc; 1 < k < M) are the numbers of clusters which have the same numbers of particles and k denotes 
the number of different cluster types. For a given k 

SI + S2 + ... + Sk = M (54) 

and 

M k 

In this representations all the integers {m^} are assumed to be different: 

1 < TOl < 7712 < ... < ruk (56) 

Due to the symmetry with respect to the momenta permutations inside the subsets of equal n's it is sufficient to 
consider the wave functions in the sectors 

qi <q2 < ... < qsi ; (57) 

qsi+i < 951+2 < ..• < qsi+s2 ; 



95i+...+sfc_i+i < qsi+...+Sk-i+2 < ... < qsi+...+Sk-1+Sk 

It can be shown (sec Appendix C) that any two wave functions 'i'q^-*(x) and '^[^Jii'x.) in which the parameters (q, n) 
and (q',n') are assumed to have the structure described above (Eqs.(l55 |) - ([571) '). are orthogonal: 

n + oo / M ^ / M 



dxi...dxAr*(*^.)(x)*J3*(x) = S{M,M')(l[S{n^,n'j)(l[{2n)5{q^-q'j) 



(58) 



where 6(n^m) is the Kronecker symbol and (5(5) is the (5-function. The above orthonormality condition defines the 
normalization constant 



C(*^)(q,n) = 



1 (n„!)2K" 



n 



^! VCrU VCrU [tea - q^ - ^-(^ ~ T + r')]^ + t\ 



1/2 



(59) 



In other words, the wave functions, Egs. ljSip or (1521) form the orthonormal set. Although, at present we are not able 
to prove that this set is complete, the suggestion of the completeness (which assumes that there are exist no other 
eigenstates besides those described above) looks quite natural. 

Finally, substituting Ea.(|Tf l) - P5)) into Eg. (1551) . for the energy spectrum one easily obtains: 



M , M 9 A/ 



= E('^^")' = -<^'& - ^E(-«-"") (60) 



^ a=l r=l ' a=l ^ a=l 



IV. REPLICA PARTITION FUNCTION 



The time dependent solution 4'(x,t) of Eg. ipsj) . satisfying the starting condition Eg. ([25)1 can be represented in 
terms of the linear combination of the eigenfunctions \I'q^n''(x), Eq. ([50l) - (|52|) : 



^^(x,^ = E E' /'^q *^'n'(x)*^^^.)*(0) exp[-i?M(q,n)i] (61) 



M=l n 
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Here the summations over are performed in terms of the parameters {si,mi}, Eas. (j54p - (j56|) : 

AI oc oo / k \ / ^ \ 

E' - E E E '^(E^«' ^(E^^-'' ^) (62) 

where for simphcity, due to the presence of the Kronecker symbols the summations over rrii and Si are extended to 
infinity. The symbol J I?q in Ea. (|6ip denotes the integration over M momenta q^, in the sectors, EQ. (j57p : the energy 
spectrum Em (q, n) is given by Eq. 

Now according Ea. (j29p for the replica partition function of the original directed polymer problem we get 



Z{N,L) = vl/(0;L) = T! /'^^ l^^^n^O)!' exp[-i?M(q, n) i] 



(63) 



M=l n 



where due to the presence of the Kronecker symbols in Eq. (|62l) the summation over M can also be extended to infinity. 
Using Ea. (|50l) . and taking into account antisymmetry with respect to the momenta permutations, one can easily prove 
that for M > 2, 

P a<b ^ 



N 



= C^^^^ (-1)1^1 n Ua(a)-g.W-»4 '° 2 (64) 

P a<b ^ 

a{a)^a{b) 

Given the antisymmetry of the product with respect to permutations of the momenta, it is sufficient to consider only 
one (trivial) permutation and multiply the result by the total number of permutations: 

M no, nfi . 

*ifj(o) = c^i^ , , • , n n n ( - '^/^ - -'+'V (^5) 

a<p r— 1 r' — 1 ^ ^ 

For M = 1, according to Eqs. (li2)) - (li5|) . 



We see that the expression |^q^n (0)|^ exp[— iJ^/ (q, n) L\ in Ea.(|551) is symmetric with respect to the permutations 
of the momenta belonging to the clusters with the same number of particles. In this case the expression for the 
partition function can be written in terms of the unconstrained integration over the momenta qa '■ 



M— 1 A;— 1 si-.-Sfc — 1 l<mi...<mi, ' ^a— 1 ^ 



oo AI oo oo r/'v^ „ n /r\ s „ „^ t\t\ / ^'i 

Z{N,L) = 

M—1 k—1 si...Sfc — 1 l<mi...<mfc 

(67) 

where 



n = {mi,..., mi, 7712,..., m2, , m^, mfe} (68) 

Eq. ([ST]) contains the summations of the quantity 

/(ni,n2,...,nM) = ( fl ^) (0)1' e-^^^-^^^ (69) 



which is the function of M integer parameters Ua- Using explicit expressions. Eg. ([BO)) . ([5^ and (I55|) one can easily 
prove that this function is fully symmetric with respect to permutations of all its M arguments. In this case 

EE E ''^■'""!:!Sr""' = M E ■- E <E".. /w (™) 

k—lsi...Sk—l l<?ni...<mfc ni— 1 riA/— 1 ^a— 1 ''^ 
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so that the summations in the expression for the replica partition function, Ea. ((67)) . can be essentially simplified: 

OO y _A/ /> + 00 ,\00 00 oo / \ 

ziN,L) = j: ^(n / ^) E E - E <E"- ^) |*(-Ho)re--(-)- (n) 

Substituting here the explicit expressions for (0), Eqs (|65|) - (|66|) . for i?M(q, n), Eqs. (|60| . as well as for the nor- 
malization constant C'^^^-' (q, n), Eg. (|59|) . one gets 

Z{N.L) = iV! Q-W'i'+^i^'-N) + 



r M oo 



M=2 ^a=ln^ = l-' " ^ 



- 9/3 - ^(J^a - - 2r + 2r') 



M ria "/3 



(72) 



The first term in the above expression is the contribution of the ground state (Af = 1), and the next terms {M > 2) are 
the contributions of the rest of the energy spectrum. Next, after a few lines of slightly cumbersome transformations 
(see Appendix D) Eg. ((72)) can be reduced to the form (cf. Eq. lfTOl) ) 

(73) 
(74) 



Z{N,L) = N\k^Q-^^^^° Z{N,L) 
where /o the linear (selfaveraging) free energy density, 

1 



and 



ZiN,X) = 



47r 



oo -j^ oo - AI 

E ^ n 



M=2 
M 

n 



-j-oo 



ni ...riM — l "-a— 1 



Pa-P0 - iX{na - np)\ 



dPa 

in 



+ 00 



M 



a=l 



a<l3 \Pa -Pl3- iX{na + np)\ 

Here instead of the system length L we have introduced a new parameter 



ML) 



1 /L 



1/3 



-(/3Vl 



1/3 



(75) 



(76) 



Next, we linearize the terms cubic in Ua in the exponentials of Eq. ()75p with the help of Airy functions, using the 
standard relation 



1 /■+°° 
exp{-X^n^) = / dy Ai(y) exp(An) 

J-oo 



(77) 



After shifting the Airy function parameters of integration ya ya + ta + Pa the expression for Z{N, A) becomes 
sufficiently compact: 



Z(iV, A) 



+ 00 



dy 

twin 



+°° dp 

-oo 47r 7o 

r M ^+oo 



dt Ai(t/ + i + p^) exp(AA^j/) + 



AI'. 

M=2 '-Q=l 



ni . . .riA/ — 1 ^a— 1 



E ^(E'^-'^in 

a<0 



oo 
M 



dpa 

47r 



+ 00 



dta Ai{ya +ta+ p^) 



Pa-Pp- iA(?- 



1 2 M 



\Pa - PP - iK'^a + np)\ 



exp(Ana?/„) 



(78) 
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Finally, after performing summations over {ni, hm} (see Appendix E) the above expression can be represented as 
an analytic function of two parameters: XN and A: 

oo 

Z{N,X) = Zi{XN) + J2 Zm{XN-X) (79) 

M=2 

where 

Zi(A^) = / dy f- / Ai(y + i+p2) exp(A7Vy) (80) 

J -00 J ^00 47r 7o 

and 

^m(AA^;A) = ^^^^ y PM(y,P) y aM(p; ^,0,x) e^-(-+-' 11 ^^(,.+.,0 _ 



with the definition 

M 



(82) 



Above we have introduced the integration operator 



/^i /•+00 ^+00 7 ^+00 
2?M(y,p) = ny '^y^'j 17 y^ dt„Ai(y„+i„+p2 



(83) 



as well as the integro-differential operator 



r M ^ ^+00 



n 



27r 9xa^ 



where it is assumed that the derivatives over {xa^} are taken at Xafi — 0. 

The crucial point is that all these factors: / T>M{y,p), J GmIp', ipifPiX), Va{'4',4>,x) as well as the last product 
in Ea. (|5T|) do not contain the replica parameter N; the latter enters the expression only in the combination AA^ 
within the exponentials of Eas. ((5n|) - (I5T]) . Thus, we have obtained the exact expression for the partition function, 
Eas. (|79l) - (|81l) . in the form of an analytic function of the replica parameter N, which until now was assumed to be an 
arbitrary integer. In the following, we will consider the analytic continuation of this function to arbitrary complex 
values of A^ and, in particular, in the limit N 0. Unfortunately this crucial step of the analytic continuation is 
ambiguous, as our partition function grows as exp(A^'^) at large A^ (it is well known that in this case there can exist 
many distribution functions which have the same moments). We will return to this problem with a further discussion 
in Section VII. 



V. THERMODYNAMIC LIMIT 



Assuming now that the parameter A^ is an arbitrary complex quantity we are going to take the thermodynamic 
limit L — T' cx) in the replica partition function, Eqs. ([75|) -(|5T |) . where, according to Eg. ([75]) . the length L enters in terms 
of the parameter X{L) = i(/3^M^L)^/'^. It is crucial that the replica parameter A^ appears in Eqs. ([79)) - (l8T]) only in the 
combination AA^. Keeping in mind further calculations of the free energy distribution function, Eqs. ([T2|) -(fT8 |) . in the 
limit A — > cx) we have to keep the value of the parameter 

s = AA^ = N (85) 

finite. In other words, in the thermodynamic limit L — > 00, the replica parameter A^ ^ 1/A ~' L^^^^ 0. 
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It turns out that in the thermodynamic hmit the expression for the rephca partition function, Eas. ([79|) - (l8T]) . 
simphfies dramatically. Indeed, since the functions 77^, Ea. (j82p . take only real values, we have 



lim 

A— >oo 



exp 


HVa 






exp 






exp 





, for ya < yi 
-1 , for ya > yi 



(86) 



Substituting this into Ea. (|5T|) we get 



M -| „ 



It can be shown (see Appendix F) that the last term in the above equation is unity 



Thus, introducing the function 



$fx) 



dp 



dt Ai(x + i 



and substituting the definition of the operator /2?j\/(y, p). Eg. (|83l) . into Eq. (f87| one gets 



Zm{s) = 



(-1) 



M-l /• + 00 



(M-1)! J_ 

Now, substituting Eqs.(IHni) and into Eq.(17ni) one finds 



dx $(x) e" 



dy^(y) 



Af-l 



lim Z{N,\) = Z{s) 

A— >oo 



The summation of this series leads to the result 



(-1)' 



M-- 



dy^{y) 



Z{s) 



dx ^(x) exp sx 



+00 



dy$(y) 



(87) 



(88) 



(89) 



(90) 



(91) 



(92) 



with the function $(x) is defined in Eq. ((89)) . Note that the prefactor A^! in the expression for the full replica 
partition function, Eq. ()73p . is irrelevant in the thermodynamic limit A 00, since for fixed parameter s — NX, 



N\ = r(l + N)n^ = ^(^ ^ ^) 



1 



A— >^oo 



and hence 



lim 

L^OG 



z{N, L) e 



fiNLfo 



Z{s) 



(93) 



(94) 



VI. FREE ENERGY DISTRIBUTION FUNCTION 



The distribution function of the free energy fluctuations can now be derived following the lines of the general 
approach discussed in the Introduction, Eas.((5|)- (I17I) . According to the definition of the replica partition function, 
Eq.®, 



Z{N,L) 



+ 00 



dF vl{f) e 



-PNF 



(95) 
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where Pl{F) is the distribution function of the total free energy of the system. Redefining the repUca partition 
function according to Ea.(|73p and introducing rescaled free energy fluctuations / according to the definition 



F = foL + -A/ 

(where /o and A are defined in Eqs. (|74|) and (|76)) ) instead of Ea. (|95|) we get 

z{N,L) = / df pune-^^'f 



where PlU) is the distribution function of the free energy fluctuations, which is related with Pl{F) via 



Taking now the limit L — >■ cx) in both sides of Ea. (P7)) . at fixed AA^ = s, we obtain 

^+00 



Z{s) 



dfp.{f)e 



where Z{s) is given in Eq dMI) . and 



P.(/) = hm PlU) 



(96) 



(97) 



(98) 



(99) 



(100) 



is the a universal thermodynamic limit distribution function of the free energy fluctuations. This function is obtained 
from the relation, Ea. ((M)) . via inverse Laplace transform 



P*if) 



ds ~ 



2m 



- z{s) e 



sf 



Substituting here the explicit expression for Z{s), Ea (j92p . we find after a trivial integrations, 

^+00 



n (/) = $(-/) exp 
where the function $(2:) is defined in Eg . (f89| . 



/ 



dy <i>iy) 



(101) 



(102) 



a, 0.1 




FIG. 1: The probability distribution function P*(/), ea. (fT02l 
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This is the central result of the paper. The shape of the distribution function P*{f) is represented in Fig.l. One 
can easily check that the function P^^if), Eq. (|102p . is positively defined and properly normalized. Indeed, the function 
^(x), Eg. (|89| ■ is positive for all values of x (—00 < a; < +00) and its asymptotics are 



>+cx)j ~ exp[~— a;' 
-00) - _ 

Thus, according to Eq. (|102p . 

/+00 r+oo J 

dfP.if) = df-exp 
-00 J-00 dj 

The asymptotics form of the left and the right tails are 



3/2 







dy $(y) 



n(/^-oo) ^ exp(--|/|3/2 



-P*(/->+oo) ~ exp -Co/' 



f3/2 



(103) 
(104) 

(105) 
(106) 



where cq — 0.11. 



VII. DISCUSSION 



Obtained result for the distribution function of the free energy fluctuations, Ea. (ll02p . is rather surprising. At 
present there arc exists an appreciable list of statistical systems for which similar distribution functions have been 
computed exactly in the thermodynamic limit. These systems are: the polynuclear growth (PNG) model [l^, the 
longest increasing subsequences (LIS) model [l3|, the longest common subsequences (LCS) [ll|, the oriented digital 
boiling model [la], the ballistic decomposition model [I4I, and finally the zero-temperature lattice version of the 
directed polymers with a specific (non-Gaussian) site-disorder distribution [l^. It is remarkable that in all these 
systems (which are not always look similar) the fluctuations of the quantities which play the role of "energy" are 
described by the same distribution function, the so-called the Tracy- Widom (TW) distribution [T2j . 

The result obtained in the present work is not the the Tracy- Widom distribution, although it is widely believed that 
the system considered here belongs to the same universality class as the models listed above. It is not clear, however, 
which way the term "universality class" should be understood. As far as only the scaling exponents are concerned, all 
these systems including the present one, indeed belong to the same universality class. On the other hand, the shapes 
of the distribution functions of the systems listed above and the result obtained for the present system are different. 
What is really surprising is that they are different not just at the level of "details" but in the qualitative properties, 
such as the asymptotic behavior of the left an right tails. The TW distribution Prwif) is strongly asymmetric in its 
asymptotics: 



Prwif ^-^) ~ exp(--|/|3/2 



Prwif ^+00) ^ exp( .f 



(107) 
(108) 



hence, the right tail decays much faster than the left one. In contrast, our distribution function is "almost symmetric" 
in its right and left tails, cf. Eq. (|105p - (|106p . and moreover, due to numerical factors its left tail decays slightly faster 
than the right one. 

One may propose three possible explanations of the discrepancy discussed above: 

(1) All the systems described by the TW distributions are essentially the zero-temperature models, while our system 
by its definition is the "high-temperature" one. Moreover, formally, our system has no zero-temperature limit at all. 
To study the limit T — >• 0, one would need to introduce a lattice or a "finite width" (5-function regularization in the 
model, cf. Eqs.([T])-([2]). In both cases, the Bethe anzats solution leading to the result, Eq. (ll02p . would no longer 
be valid. On the other hand, from a general physical point of view the conclusion that the thermodynamic limit of 
directed polymers at T = and in the limit T — look different would be quite surprising. In particular, the solution 
in the thermodynamic limit of the present system confined to a cylinder geometry reveals no such difference [l9j . 
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(2) Unlike the present system, the disorder in all the systems described by the TW distribution is essentially non- 
Gaussian. Formally, any deviation from the Gaussian statistics of the disorder would again ruin our Bethe anzats 
solution, and at present it is not clear to what extend this solution is "stable" with respect to non-Gaussian deviations. 

(3) Much more likely appears the technical problem given by the third option. According to Eq. ([M)) the distribution 
function P»(/) is defined by the replica partition function Z{s) at finite (both positive and negative) values of the 
parameter s cx L^^^N, which means that in the thermodynamic limit L — >■ oo the function P*{f) is defined by the 
replica partition function with the "number of replicas" iV — >■ (both positive and negative). In the present work 
we have computed the replica partition function for arbitrary positive integer N with the aim to perform an analytic 
continuation for the values N 0. Unfortunately, the analytic continuation from integers to arbitrary (real or 
complex) N is unambiguous only if the corresponding function growth at infinity not faster than exp(A^). In our case 
the partition function growth as exp(A^^), and hence the knowledge of this function at arbitrary integer N does not 
guarantee an unambiguously reconstruction in the region |iV| ^ 1 [11| . The classical example of such type of situation 
is well known in the theory of the mean-field spin-glasses, such as the Sherrington and Kirkpatrick (SK) model [2^ 
and the Random Energy Model (REM) For both models one can relatively easy compute the replica partition 
functions for arbitrary positive integer number of replicas N, and in both cases the replica partition function growths 
as exp(iV^) at large N. Further "direct" analytic continuation of these solutions to the region < < 1 yields 
nothing else but the replica-symmetric solutions, which at first sight look sufficiently reasonable (at least in the SK 
model), but more detailed investigation reveals that they are unphysical. As we know, the solution which is believed 
to be valid in the region < A < 1 reveals the Parisi replica symmetry breaking (RSB) structure (one-step RSB in 
the case of REM), and it is derived in terms of a heuristic procedure (directly in the interval < A < 1) and not as 
a proper analytic continuation from integer to non integer values of A". Moreover, in the case of REM there is a kind 
of the phase transition at A^ = 0, which means that at negative A^ the replica partition function should be computed 
separately [31]. 

Although up to the present moment, we have not uncovered any unphysical properties in obtained probability 
function P*{f), the above arguments indicate that the present solution, Eq. (|102p . could as well be a kind of a distant 
analog of the "replica symmetric approximation", while the derivation of the "true" solution would require more 
sophisticated ideas. 
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Appendix A 

Wave functions of quantum bosons with repulsive interactions 



Explicitly the eigenstate equation (I5T|) reeds 



N 1 ^ 



2/3 



(A.l) 



Due to the symmetry of the wave function with respect to permutations of its arguments it is sufficient to consider it 
in the sector 



Xi < X2 < ■•■ < Xn 
as well as at its boundary. Inside this sector the wave function 5'(x) satisfy the equation 



(A.2) 



a=l 



(A.3) 



which describes N free particles, and its generic solution is the linear combination of N plane waves characterized by 
TV momenta {gi, g2, 9^} = ^.■ Integrating Eq. (jA.l| over the variable (x^+i — Xi) in a small interval around zero, 
l^i+i — Xi\ < e — ^ 0, and assuming that the other variables {xj\ (with ^ i, i + 1) belong to the sector, Ea. (|A.2[ ). one 
easily finds that the wave function ^'(x) must satisfy the following boundary conditions: 



-5,, + k)^'(x) 



= 



(A.4) 



where k = P^u. Functions satisfying both Eq. f|A.3P and the boundary conditions Eq. (jA3| can be written in the 
form 



^q,...q^{xi,...,XN) = ^'^^\^) = C[\[[d^^-d^,+K] \ det[exp(^qx)] 



N 



a<b 



(A.5) 



where C is the normalization constant to be defined later, and the symbol exp(iqx) denotes the N x N matrix with 
the elements exp(iqaXb) {a,b = I, N). First of aU, it is evident that being the linear combination of the plane waves, 
the above wave function satisfy Eq. (jA.3p . To demonstrate which way this function satisfy the boundary conditions, 
Ea. (IA.4l) . let us check it, as an example, for the case i = 1. According to Ea. (IA.5[) . the wave function ^q^''(x) can be 
represented in the form 



*W(x) = -(5.^-5.,-«)^W(x) 



(A.6) 



where 



JV 



a— 3 



*(j^)(x) = C m - d,, + At] -d,,+K\)( II - d,, + k] det[exp(zqx)] (A.7) 



N 



3<a<b 



One can easily see that this function is antisymmetric with respect to the permutation of xi and X2. Substituting 
Ea. (IA.6l) into Eq. (|A.4p (with i = 1) we get 



{dx2 - - 



*W(x) 







(A.8) 



Given the antisymmetry of the l.h.s expression with respect to the permutation of xi and X2 the above condition is 
indeed satisfied at boundary Xi = X2- 
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Since the eigenfunction ^'q^''(x) satisiying Eg. ljA.ip must be symmetric with respect to permutations of its argu- 
ments, the function, Eg. dA.SI) . can be easily continued beyond the sector, Ea. (IA.2l) . to the entire space of variables 
{xi,X2, ■.■,xn} e Rn, 



N 



= C[Y[ - a,,) + iKSgn{xa - Xb)] det[exp(iqx)] (A.9) 



a<b 



where, by definition, the differential operators dx^ act only on the exponential terms and not on the sgn(a;) functions, 
and for further convenience we have redefined i^(^~^y^C — > C. Explicitly the determinant in the above equation is 



N 



det[exp(iqx)] = ^(-l)'^' exp [z ^ (jp„ Xa] (A.IO) 



a=l 



where the summation goes over the permutations P oi N momenta {qi,q2, ...,qN} over N particles {xi,X2, ...^xn}, 
and [P] denotes the parity of the permutation. In this way the eigenfunction, Eq. (|A.9|) . can be represented as follows 



N s N 



= C^(-l)[^l il[[-i{dx^-dx,)+insgn{xa-xb)]] exp[iY,qpa^a] (A.ll) 

P ^a<b ^ a=l 

Taking the derivatives, we obtain 

, N . N 

vI/W(x) = C^(-l)[^l (l[[qp^-qp,+insgn{xa-Xb)]] exp[i^gp„x,] (A.12) 

P ^a<b ' a=l 

It is evident from these representations that the eigenfunctions (x) are antisymmetric with respect to permuta- 
tions of the momenta qi, ...,qN- 

Substituting the expression for the eigenfunctions, Eq. ljA.sp (which is valid in the sector, Eq. (|A.2p ). into Eq. (|A.3P 
for the energy spectrum we find 

1 ^ 

^ a=l 

Now one can easily prove that the above eigenfunctions with different momenta are orthogonal to each other. Let 
us consider two wave functions ^q^-'(x) and ^'^^•'(x) where it is assumed that 

qi < q2 < ... < qw (A.14) 
q[ < q'2 < ... < q'^ 

Using the representation, Ea. (|A.lip . for the overlap of these two function we get 

/+00 
d-x<-'*(xX)(x) 
00 

.+00 ^ / JV ^ TV 



K 

Integrating by parts we obtain 



|q2^(_^)[P]+[P'] I d''J(\{[^{^x^^^x,)-^^sgn{xa~Xb)\] exp ^ x,] j X 

P,P' ^ ^a<b ^ a=\ > 

1 ( II -dxi)^ msgn(a;a - Xf,)] j cxp[i^ ] (A.15) 



/+00 ^' 
d^xexpH^,;,x.] X (A.16) 

t-.t- a=l 
,N s AT . 

X \\\\-i{dx^- - iK.%gx\(xa- Xb)\ - 9^;,) + m sgn(xa - Xf,)] j exp[i ^ ^p^Xq] \ 

^a<b ' 0=1 ^ 
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P,P' a=l ^a<b 



TV 



Taking the derivatives and performing the integrations we find 



(A.17) 



<"^*(x)<)(x) . |C|^^(-l)[-l+[-'l 

p,P' 

= |C|2^(-l)[^I+[^'l 



p.p' 



a<(, ^ a=l 

TV \ r ^ 



(A.18) 



Taking into account the constraint, Ea. (|A.14[) . one can easily note that the only the terms which survive in the above 
summation over the permutations arc P = P', all contributing equal value. Thus, we finally get 



N 



a<b 



N 



l-a=l 



With the normalization constant 



\C\ = 



(A.19) 



(A.20) 



we conclude that the set of the eigenfunctions, Ea. (|A.ll[) or (jA.12p . are orthonormal. The proof of completeness 
of this set is given in Ref. [13]. It should be noted that the above wave functions present the orthonormal set of 
eigenfunctions of the problem, Ea. (IA.l[ ). for any sign of the interactions k, e.i. both for the repulsive, k < 0, and for 
the attractive, k > 0, cases. However, only in the case of repulsion this set is complete, while in the case of attractive 
interactions, k > 0, in addition to the solutions, Eq. ljA.lip . which describe the continuous free particles spectrum, 
one finds the whole family of discrete bound eigenstates (which do not exist in the case of repulsion) (see Appendices 
B and C). 



Appendix B 

Ground state of quantum bosons with attractive interactions 



The simplest example of the bound state is the one in which all N particles are bound into one finite size "cluster": 

N , N 



*W(x) = c exp 



iq 



-K ^ \Xa 



Xb\ 



,,6=1 



(B.l) 



where C is the normalization constant (see below) and q is the (continuous) momentum of free center of mass motion. 
Substituting this function in Ea. (|A.ip . one can easily check that this is indeed the eigenfunction with the energy 
spectrum given by the relation 



E = -- 



1 ^ r 1 ^ 



b=l 



(B.2) 
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where it is assumed (by definition) that sgn(O) = 0. Since the result of the above summations does not depend on the 
mutual particles positions, for simplicity we can order them according to Eq. (jA.2p . Then, using well known relations 



N 



6=1 



N 

a=l 

JV 

a=l 



-(iV + 1 - 2a) 
N{N + 1) 
(A^ + l)(2iV+l) 



for the energy spectrum, Eq. (|B.2|) . we get 

E = 



The normalization constant C is defined by the orthonormahty condition 



/-i-oo 
dxi...dxN ^[1^ (x)*^^)(x) = i2n)Siq-q') 
-OQ 

Substituting here Eq. ljB.ip we get 

^'^}^*(x)*^^^(x) = |Cp / dxi...dxN exp i{q - q')^Xa - -^^ E ^ ^ 



(B.3) 
(B.4) 
(B.5) 

(B.6) 
(B.7) 



a=l 



a,fc=l 



r+oo n+oo 


r+oo 




/ dxi / dx2--- 


. dxN exp 


i{q 


J —OQ J Xi 


J Xn-1 





N 



N 



iiq - q')J2^a + nj^i^ + - 2a)a;a 



(B., 



where for the ordering, Eq. (jA.2[) . we have used the relation 



1 ^ 

- \xa-Xb\ = -Y,iN +l-2a)xa 



(B.9) 



a,b—l a—1 

Integrating first over xn, then over xn~i, and proceeding until xi, we find 



*(^>*(x)vi/W(x) = \cfm n 



1 



;Li r[{N ~r)K~i{q- q' 

N-l 

^ {27r)6{N{q-q')) 



--] I dxiex.]i[iN{q- q')xi] 



1 



II I 11 r(N-r)K 
Nk 

According to Eg. (IB. 71) this defines the normalization constant 



\C\ 



kN 



(B.IO) 



(B.ll) 



Note that the eigenstate described by the considered wave function, Eq. (jB.l[) . exists only in the case of attraction, 
K > 0, otherwise this function is divergent at infinity and consequently it is not normalizable. 

It should be noted that the wave function, Eq. (|B.ip . can also be derived from the general eigenfunctions structure, 
Ea. (IA.12[) . by introducing (discrete) imaginary parts for the momenta qa- We assume again that the position of 
particles are ordered according to Eq. (|A.2p . and define the particles' momenta according to the rule 



qa = q~-KiN+l-2a) 



(B.12) 
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Substituting this into Ea. (|A.12[ ) we get 

Vf^lxi < X2 < ... < X^) OC ^(-1)[^1 f n 

P ^a<b 

r N ^ N 

xexp iq^Xa + -^{N + l-2Pa)Xa 



N 



{q - '-n{N + 1 - 2Pa)) - {q - '-k{N + 1 - 2Pb)) - ik) 



, N N N 



a<b 



a=l 



(B.13) 



(B.14) 



Here one can easily note that due to the presence of the product Y\!^^i,[Pb — Pa + ^] in the summation over permutations 
only the trivial one, Pa = a, gives non-zero contribution (if we permute any two numbers in the sequence 1, 2, N 
then we can always find two numbers a < b, such that Pf, = — 1). Thus 



'^'^J'^ixi < X2 < ... < xn) oc exp 



N ^ N 

iqY.Xa+'^Y.'^N + l- 2a)xa 



a=l 



a=l 



(B.15) 



Taking into account the relation, Eq. (jB.9|) . we recover the function, Eq. (jB.l[) . which is symmetric with respect to its 
A'' arguments and therefore can be extended beyond the sector, Eq. (jA.2p . for arbitrary particles positions. Finally, 
substituting the momenta, Ea. (|B.12p . into the general expression for the energy spectrum, Ea. (|A.13p . we get 



N 



(B.16) 



Performing here simple summations (using Eqs. (|B.4p . (IB.Sp ) one recovers Eq. (|B.6p . 



Appendix C 

Wave functions of quantum bosons with attractive interactions 



1. Eigenfunctions 



The general expression for the eigenfunctions both for the case of repulsion and for the case of attraction is given 
in Eas. ([55| or ([57)1 - (1551) . A generic eigenfunction is characterized by N momenta parameters {qa} {a = 1,2,...N) 
which in the case of attractive interactions may have imaginary parts. It is convenient to group these parameters into 
M {1 < M <N) "vector" momenta, 

= - ^ K (n„ + 1 - 2r) (C.l) 

where qa {a = 1,2, M) are the continuous (real) parameters, and the (discrete) imaginary components of each 
"vector" are labeled by an index r = 1, 2, n„. With the given total number of particles equal to N, the integers 
have to satisfy the constraint 

M 

J2 na ^ N (C.2) 

In other words, a generic eigenstate is characterized by the discrete number M of complex "vector" momenta, by the 
set of M integer parameters {rii, n2, tim} = n (which are the numbers of imaginary components of each "vector") 
and by the set of M real continuous momenta {qi, q2, qn} = Q- 
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To understand the structure of the determinant of the N x N matrix exp(i qa xi,), which defines the wave functions, 
Eas. ([57| - p8p . the N momenta Qa = q" can be ordered as follows: 

{qa} - {ql ql ; ql ql ql,; , q^\ (C.3) 

By definition, 

N 

det[ex^^x)] = XI exp[i^qp^xa] (C.4) 

P a=l 

where the summation goes over the permutations of N momenta {qa}, Ea. (|C.3p . over N particles {xi, X2, xat}, and 
[P] denotes the parity of the permutation. For a given permutation P a particle number a is attributed a momentum 
component q^^^y The particles getting the momenta with the same a (having the same real part qa) will be called 
belonging to a cluster f^Q,. For a given permutation P the particles belonging to the same cluster are numbered by 
the "internal" index r ~ 1, ...,nQ. Thus, according to Eq. ((38|) . 



N 



H'n^x) = C^^)X(-1)[^1 ^[-^(a.^-a,J+^.sgn(x<,-.,)] exp 



a<b 



N 



(C.5) 



where Cq^n'' is the normalization constant to be defined later. Substituting here Ea. (|C.l[) and taking derivatives we 
get 

N r 



a<b 



qa{a) - qa{b) ~ ^n[ ^ r{a) + r[b) - sgn{Xa ~ Xb)\ 



X exp 



AT 



a=l 



N 



i^qa{a)Xa + -^^{^aia) + 1 - 2r{a))xa 



(C.6) 



The pre-exponential product in the above equation contains two types of term: the pairs of points (o, b) which belong 
to different clusters (a(a) ^ and pairs of points which belong to the same cluster (a(a) = a{h)). In the last 

case, the product lia over the pairs of points which belong to a cluster Via reduces to 

Ha (X n [r{b)-r{a)-sgTi{xa-Xb)]] (C.7) 

a<b&Q,a 

As for the ground state wave function Eq. (jB.14p - (|B.15[) . one can easily note that due to the presence of this product 
in the summations over Uq! "internal" (inside the cluster Via) permutations r(a) only one permutation gives non- 
zero contribution. To prove this statement, we note that the wave function 5'q^''(x) is symmetric with respect to 
permutations of its N arguments {xq}; it is then sufficient to consider the case where the positions of the particles 
are ordered, a:i < X2 < • • • < xn- In particular, the particles {xa^} (^ = 1; 2, . . . , Ua) belonging to the same cluster 
Q.a are also ordered Xa^ < Xa2 < ■ ■ ■ < Xa„^ ■ In this case 



n„ cx J][r(0 -r(fc) + l] 



(C.8) 



k<l 



Now it is evident that the above product is non-zero only for the trivial permutation, r{k) — k (since if we permute 
any two numbers in the sequence 1, 2, ...Tna, we can always find two numbers k < I, such that r{l) = r{k) — 1). In 
this case 



(C.9) 



k<l 



Including the values of all these "internal" products, Eq. (|C.9|) . into the redefined normalization constant C^n , the 
wave function, Eq. (|C.6|) (with xi < X2 < ■ ■ ■ < xn), reeds 



N 



a<b 
i(a)#a(b) 



qa{a) - qa(b) " «k[ r{a) + r{b) + 1\ 



X exp 



AT N 



(C.IO) 
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where the product now goes only over the pairs of particles belonging to different clusters, and the symbol X]p' means 
that the summation goes only over the permutations P in which the "internal" indices r(a) are ordered inside each 
cluster. 

Note that although the positions of particles belonging to the same cluster are ordered, the mutual positions of 
particles belonging to different clusters could be arbitrary, so that geometrically the clusters are free to "penetrate" 
each other. In other words, the name "cluster" does to assume geometrically compact particles positions. 

Now using the symmetry of the wave function 5'q*n (x) with respect to the permutations of its arguments the 
expression in Eq. ljC.lOp can be easily continued beyond the the sector Xi < X2 < ■■■ < xn for the entire coordinate 
space i?jv. Using the relations 



^ (n„ + 1 - 2r(a)).Ta = + 1 - 2fc)a;„, = XI 1^'^'= 



(C.ll) 



k,l = l 



and 



(na + 1 - 2fc) = - Xsgn(a;aj, ~ 



(C.12) 



1=1 



(where Xa^ < < ■■■ < Xa„^), the wave function ^q^n (x), Ea. (|C.10[) . with arbitrary particles positions reeds 



N 



a<b 
a{a)^a{b) 



qa{a) - qa{b) + Y X! ^gn(^a ~ ^c) ~ y ^ Sgn{xb - X^) + IK Sgn{Xa ~ Xb) 



(i>) 



X exp 



r M M 



^ a— 1 aeOc 



a— 1 a.a'eSla 



(C.13) 



Here the summation goes only over the permutations P of the momenta, Eq. (jC.3l) . in which "internal" indices r{a) 
in the clusters are ordered according to the actual spatial ordering of the particles belonging to these clusters (i.e. 
r(a) increases from the smallest Xa in the cluster to the largest one). This wave function can also be re- written in the 
more compact form: 



N 



a{a)^a{b) 



-i{dx^-dxt) +iKSgIl{Xa-Xb)] 



exp 



At 



M 



a— 1 aG^^c^ 



Q — 1 a^a'^^la 



E E 4 E E 

(C.14) 



2. Orthogonality 

We define the overlap of two wave functions characterized by the two sets of parameters, (M, n, q) and (M', n', q') 



as 



qLT'M) - / '^^x^ffj (xX^)(x) 

^ — oo 

Substituting here Ea. (|C.14l) we get 

p pi J~oo 



N 

n 

a<b 
a' {a)^a' {b) 



N 

n 

a<b 
(y{a)^a{b) 



i{dx^ - dx^) - iKSgn(xa - Xb) 



'i{dx^ - dxb) + iKSgn{xa ~ Xb) 



exp 



M' "o 



E 9« E - 1 E E 1^'^ 



exp 



Q — 1 aGO' 



M 



a— 1 a,b£fl' 



M 



^E^" E ^--jE E 



a— 1 a^Qck 



a— 1 a,bGCla 



(C.15) 



(C.16) 



2^6 



\Xa - Xb\ 
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where {ila} and {^2^} denote the clusters of the permutations P and P' correspondingly. Integrating by parts we 
obtain 



p pr J —CO 



a x exp 



N 

n 

a<b 



iKSgll{Xa - Xb) 



M' 



E 9" E ^'^ - I E E 1^'^ 



Q— 1 aEO 

n 



a=l a,befl' 



'iipx^ ~ dx^) + iKSgn(a;a - Xb) 



X exp 



M M 

: E 9" E - I E H 1^- 



a;6 



a— 1 a.^Gfio 



(C.17) 



First, let us consider the case when the integer parameters of the two functions coincide, M = M' and n — n', and 
for the moment let us suppose that all these integer parameters {na\ are different, 1 < rii < n2 < ... < um- Then, 
in the summations over the permutations in Eq. (|C.17[) . we find two types of terms: 

(A) the "diagonal" ones in which the two permutations coincide, P = P' ; 

(B) the "off-diagonal" ones in which the two permutations are different, P ^ P' . 
The contribution of the "diagonal" ones reeds 



Q 



d^xexp-^E9:«E-'^-iE E K 



Xb\ 



N 

n 



a— 1 a£Vt 
r M 

exp 



M M 

E E - 1 E E 



a— 1 a.b^VLa 



(C.18) 



It is evident that all permutations a (a) in the above equation give the same contribution and therefore it is sufficient 
to consider only the contribution of the " trivial " permutation which is represented by line in Eq. (jC.Sp . The cluster 
ordering given by this permutation we denote by Q;o(a). For this particular configuration of clusters we can redefine 
the particles numbering, so that instead of a "plane" index a — 1,2,...,N the particles would be counted by two 
indices {a,r): {xa} {^r} — M) {r — indicating to which cluster a a given particle belongs and 

what is its "internal" cluster number r. Due to the symmetry of the integrated expression in Eq. (|C.18|) with respect 
to the permutations of the particles inside the clusters, we can introduce the "internal" particles ordering for every 
cluster: xf < X2 < ... < x"^. In this way, using the relation, Eg. dC.llj) . we get 



Nl 



M 

n 



+00 



+ c>o 



dx^ 



-^00 



X exp 



X 



M rio 



E^^E^" + 2 E Ek + 1 - 2^)< 

a— 1 r— 1 a— 1 r— 1 

M ria p 

nnn -(5.^-9../+. 

z E 9" E < + i E E(^^" + 1 - 



^a<fl r=l r' = l 

M n. 

X exp 



Q — 1 r— 1 



(C.19) 



where the factor N\/ni\...nM^- is the total number of permutations of M clusters over N particles. Taking the 
derivatives and reorganizing the terms we obtain 



[qa - - —{ua - np - 2r + 2r')] 



2 ? 



^a</3 r—1 r' — l 

( r-\-oo r-{-oo r-\-oo | 

^riW ^""i dx^.... da;"^ e'('?°-«")^"=i<+''^"=i^"°+^-^'')< WC.20) 

Q = l [ "'-00 -^X? •^"^"a-l J 
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Simple integrations over yields (cf. Eqs. (IB.8l) - (|B.lip ') 



a<f) r=l r' = l 



AI ria 



[<?a - 90 - Y (j^a - n-^ - 2r + 2r')] ^ + 



n 



(C.21) 



Now let us prove that the "off-diagonal" terms of Ea. (|C.17p . in which the permutations P and P' are different, 
give no contribution. Here we can also chose one of the permutations, say the permutation P, to be the "trivial" one 
represented by line in Eq. (jC.3p with the cluster ordering denoted by ao(a). Given the symmetry of the wave functions 
it will be sufficient to consider the contribution of the sector xi < X2 < ■■■ < Xn^- According to Eq. (jC.17l) . we get 



QiV'^'^'M) oc Y^{-ir' / d^x exp ^ ^ \x. 

P' J xi<...<x„^ L a=l aefl'^ a = la,ben'^ 



Xb\ 



N 

n 

a<b 



-i{9x^ - dxb) + iKSgn{xa - Xb) 



N 

n 

a<b 



-i{dxa - dx^) - iKSgn{xa - Xb) 



X exp 



M M 

' X! 9" X! + 9 X! X! + ^ 2r{a))xa 



Q— 1 a^Q.° 



Q — 1 a^n° 



(C.22) 



Here the symbols {f2° } denote the clusters of the trivial permutation ao(a). Since P' ^ P, some of the clusters ft'^ 
must be different from il° . As an illustration, let us consider a particular case of N = 10, with three clusters ni = 5 
(denoted by the symbol "O") , ^2 — 2 (denoted by the symbol " x ") and ria = 3 (denoted by the symbol 'A"): 



particle number a 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


permutation ao{a) 

















X 


X 


A 


A 


A 


permutation a' (a) 











A 





X 


X 





A 


A 



Here in the permutation a' (a) the particle a = 4 belong to the cluster a = 3 (and not to the cluster a = 1 as in 
the permutation ao(a)), and the particle a = 8 belong to the cluster a = 1 (and not to the cluster a = 3 as in the 
permutation ao{a)). Now let us look carefully at the structure of the products in Ea. (IC.22p . Unlike the first product, 
which contains no "internal" products among particles belonging to the cluster U^, the second product does. Besides, 
the signs of the differential operators (^dx^ — dx^) in the second product is opposite to the "normal" ones in the first 
product (cf. Eqs. (|C.7p - (IC.9p V It is these two factors (the presence of the "internal" products and the "wrong" signs 
of the differential operators) which makes the "off-diagonal" contributions, Ea. (|C.22p . to be zero. Indeed, in the above 
example, the second product contains the term 



n4,5 = 



exp 



aef2° 



a— 1 a(EO° 



(C.23) 



(we remind that the particles in the clusters 17° are ordered, and in particular x^ < x<^) 
get 



Taking the derivatives, we 



H' 



4,5 



iqi + ^ ("1 + 1 - 2r-(4)) - - I (m + 1 - 2r(5)) 



[r(4) - r(5) + l] = 



(C.24) 



since in the first cluster r(a) = a. 

One can easily understand that the above example reflect the general situation. Since all the cluster sizes Ua are 
supposed to be different, whatever the permutation a'{a) is, we can always find a cluster 12° such that some of its 
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particles belong to the same cluster number a in the permutation a' (a) while the others do not. 
consider the contribution of the product of two neighboring number points 



Then one has to 



/c,/c+l 



exp 



M 



9a J2 



M 
K \ ^ 



^ (n-a + 1 - 2r{a))xa 



a—1 aGO° 



(C.25) 



where in the permutation a' (a) the particle k belong to the cluster number a and the particle (fc + 1) belong to some 
other cluster. Taking the derivatives one gets 



fc,fc+i 



cx [r(fc)-r(fc + l) + l] = 



(C.26) 



as r(a) is the "internal" particle number in the cluster r2°, where r{k + 1) = r{k) + 1 (cf. Eqs. (|C.7|) - (|C.9j) ). 

Thus, the only non-zero contribution to the overlap, Eq. (jC.15|) . of two wave function ^'q*^^(x) and (x) (having 
the same number of clusters M and characterized by the same set of the integer parameters 1 < rii < n2 < ... < um) 
comes from the "diagonal" terms, Ea. (|C.21l) : 



M 



M Tta 



Q = l 



1/3 r 



x'nnn 

^Q</3 r=l r' = l 



[qa - g/j - Y ("a - "-/3 - 2r + 2r')] ^ + 



M 

n 



(2^)<5(q„ - q'J 



(C.27) 



The situation when there are clusters which have the same numbers of particles Ua is somewhat more complicated. 



Let us consider the overlap between two wave function ^q*^2(x) and ^'q.n''(x) (which, as before have the same M 



and n) such that in the set of M integers ni, n2, ?ia/ there are two n^'s which are equal, say 



(where 



Oil 7^ 012)- In the eigenstate (q',n) these two clusters have the center of mass momenta q'^^ and g^^, and in the the 
eigenstate (q, n) they have the momenta qa-^ and qa^ correspondingly. According to the above discussion, the non-zero 
contributions in the summation over the cluster permutations a (a) and a' {a) in Eq. (|C.17P appears only if the clusters 
{fJct} of the permutation a{a) totally coincide with the clusters {^2^} of the permutation a' {a). In the case when all Ua 
are different this is possible only if the permutation a{a) coincides with the permutation a' [a). In contrast to that, in 



the case when we have Ua-^ 



there are two non-zero options. The first one, as before, is given by the "diagonal'' 



terms with a{a) = a' {a) (so that the clusters {^0} and {51^} are just the same), and this contribution is proportional 
to 5{qa^ — ) 5{qa2 — Qui)- second ("off-diagonal") contribution is given by such permutation a' (a) in which the 
cluster (of the permutation a'{a)) coincide with the cluster fi^j (of the permutation a{a)) and the cluster f2^^ 
(of the permutation a'{a)) coincide with the cluster il^^ (of the permutation a{a)) while the rest of the clusters of 
these two permutations are the same, fl^ = fla (a ^ 01,0:2). Correspondingly, this last contribution is proportional 
to 6{qa^ ~ q'a^)5{qa2 — ^qj) (— l)"°i. In fact this situation with two equivalent contributions is the consequence of 

the symmetry of the wave function 'I'q^^,^(x): the permutation of two momenta qa^ and qa^ belonging to the clusters 
which have the same numbers of particles, ria^ = tIq^ produces just the factor (— (see discussion below Ea. ([52|) '). 
Therefore considering the clusters with equal numbers of particles as equivalent and restricting analysis to the sectors 
?ai < <la2'i I'ai < 9q2 ^^'^ thsit the sccond contribution, 5{qai — <l'a2)^{loL2 ~ Qai) identically equal to zero, thus 
returning to the above result Eq. (|C.27[) . 

In a generic case the M-component vector n can be represented in the form 



n = { mi, . . .,TOi , m2, . . .,m2 , , mfc, . . .,mfc } 

Si S2 Sfc 

where si -I- S2 + ■•■ + Sfc — M and k integers {mi} (1 < k < M) are all supposed to be different: 

1 < TOl < TO2 < ■•■ < 



(C.28) 



(C.29) 



Due to the symmetry with respect to the momenta permutations inside the subsets of equal n's it is sufficient to 
consider the wave functions in the sectors 



gi < 92 < •■• < qsi ; 

Qsi+i < <Zsi+2 < ••■ < qsi+s2 ; 



(C.30) 



9si-t-...+Sfc_i + l < qsi+...+Sk-i+2 < ■■■ < qsi+...+Sk-1+Sk 
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In this representation we again recover the above result Eq. (jC.27|) 

Finally, let us consider the overlap of two eigenstates described by two different sets of integer parameters, n' ^ n. 
In fact this situation is quite simple because if the clusters of the two states are different from each other, it means 
that in the summation over the pairs of permutations P and P' in Ea. (IC.17p there exist no two permutations for 
which these two sets of clusters {Q,a} and {il^} would coincide. Which, according to the above analysis, means that 
this expression is equal to zero. Note that the condition M' ^ M automatically implies that n' ^ n. 

Thus we have proved that 



M 



:m n 



a=l 



(r7,Q!)2K»» 



M Uc, rifj . 



nnn 

_ ^a</3r— Ir' — 1 



[qa - qp^Y -np-2r + 2r')] ^ + 



(C.31) 



where the integer parameters {ua} and {n^} are assumed to have the generic structure represented in Eas. (jC.28[) - 
(jC.29|) . and the momenta {qa} and {q'^} of the clusters with equal numbers of particles are restricted in the sectors, 
Ea. (IC.30[) . Finally, according to Ea. (jC.31|) . the orthonormality condition defines the normalization constant as 
it is given in Eq. ([5^ . 



Appendix D 



In this Appendix we simplify the last term in the r.h.s of Ea. (|72l) . and prove that 

j-j ^ -Q -Q \qa -qp- —{ng - n/j - 2r + 2r )| _ |ga - qp ~ —{ng - np)] 
r=ir'=i[qa-qi3~ f{na~np-2r + 2r')] +k2 \qa-- qp - ^{na+npY 



We rewrite 11 as 

n":i U7=i Ife - f rv) - {qp - fnp) + in{r - r')\ 



n 



nr=i [(9a - + inr) - {qp - fnp + inr') - in] [{q^ - fn^ + inr) - {qp - fnp + inr') + in] 

(D.2) 

Redefining the indices r and r' of the product in the left brackets [...], of the denominator 

r Ua + l-r (D.3) 
r' -)■ np + l-r' 

we find 

n ^ n":i n"ii |(ga - fng) - {qp - f n^) + ^n[r - r'f 

n":i KU [(9a + f "a) - ilP + T^P) - *«(^ - ^' + 1)] [(9a " f " (^Z? " f "Z^) + - r' + I)] 



KU I (9a - f na) - {qp - f n^) + *«(r - r') 



2 



nr=i n"ii I (9a - f r^a) - (g^ - f n^) + *«(r - r' + 1)|' 
Now, shifting the product over r in the denominator by 1 we obtain 

n:=i KU |(9a - f n.) - (g^ - f n^) + ZK(r - r'f n^ti |(9a - f - [qp f n^) + z«:(l - /)|' 



n 



n"=r n"ii I (9a - f na) - - f n^) + »^(r - r')] n"ti | (9a - f - {qp - fnp) + + 1 - 

(D.6) 
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Redefining r' in the product in the denominator, 

r' ufi + l-r' (D.7) 

we obtain 

jj ^ n"ii - f n.) - (qp - f n/^) - ^<r' - 1)|' ^ KLi \ (go - f ^g) " ilp ~ fnp) " - 1)1' ^ 
Finahy, shifting the product in the numerator by (—1), we get 

n"!! - ^^a) - {qp - fnp) - inr'l^ |(g„ - fn^) - {qp - ^np) - iKfip]^ |gf„ - qp - ^{ua + np)^ 

(D.9) 

which proves Ea. (|D.l[) . 



Next, we substitute Eq. (ID.1[) into the expression Eq. ([7^ for the partition function Z{N,L) and represent it in 
the form Eq. ((73)) . where 



Z{N.L) 



+ 00 



dq 1 ^^ULa^+iili^MS 
27r kN 



1 



+ Em ^ n 



M 



+ 00 



dqa 1 



Af=2 " ' ni...riA/ = l'-ct=l °° 

M 



2tT KUa 



M 



n 



Redefining all the momenta 



and introducing a new parameter 



we rewrite: 



qa 



13k 



1/3 



2\ 1/3 



(D.IO) 



(D.ll) 



(D.12) 



2/3'^° 



(D.13) 



qa~ qp~ ^{ng ~ np) _ pp- X{na - np) 

qa- qp~ fina + np) Pa ~ pp - X{na + np) 



and obtain the integrals: 
27r Kria 



/ 



+00 



+00 



[...] . / 1^^/ dte--- 



+ 00 



dte 



— An„ i 



[...] 



(D.14) 



(D.15) 



Substituting the transformations, Ea. (ID.13P - (|D.15l) . into Eq. (|D.10|) . we get Eq. ((75)) . 
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To perform the summation over ni,...,nM in Ea. (|78p . let us modify the pre-exponcntial factor: 



M 



F(p;n) ^ n 



\Pa -Pl3 - i\{na - 



M 

n 



c<p \Pa - PP - iK^a + np)\ ^ -Pp\+ iK^c + np)) ~Pp\~ iA(n„ + np)) 

We introduce the auxihary fields Xaf}, 0q0, and ipap to generate the numerator 

d 



{\Pa -Pp\ + i\{na - np)) 
and the denominators 



dXap 



exp[(|pQ - P/3| + iX{na - ni3))xap] 



X»/3=0 



-y/ (ba - Ppl + iX{na + np)) 
1 

\J {\Pa - Pi3\ - iX{na + np)) 
EH-dEll we I 

dTpapd(j>ap d 



exp 



exp 



Combining together Eas. (jE.2l) - (jE.4p we rewrite the prefactor as 
where 

M 



;{\Pa - Pl3 \ + iXina + np))(j)lp 



-T^iW-Ppl - ^Hria + np))iplp 



M 



{Xo/3}=0 



Introducing now the integro-differential operator 

M I. I.+00 
^a^P 

as well as the integration operator 



d')papd(j)ap d 

27r 9xq^ 



(E.l) 

(E.2) 

(E.3) 
(E.4) 

(E.5) 
(E.6) 
(E.7) 
(E.8) 



the partition function, Eq.(|78p. can be represented as follows 
Z{N,\) = / exp(A7Vy) + 



Af=2 ■ ni...nM = l ^Q = l ^ a=l 



(E.9) 



Next we prove that 

oo 

S{N,M) ^ E «i 



ni...nAf = l 



a 



"1 ^"2 „n3 
"2 "'3 



as 



Af 



<r ^ E"- ^ 



as 



(ai - 02) (ai - as) (ai - aaj- ) 

AT Qi Q2 QAf 

^ (as - ai) (as - 02) "' (03 - aM) 



(02 - ai) (02 - 03) (02 - flAf ) 

AT fll ^2 OAf-l 



*Af 



(flAf - fll) (flAf - 02) (flAZ-aAf-l) 



(E.IO) 
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To unbound the summations over ni, 712, riM in the above series, let us introduce the integral representation of the 
Kronecker symbol: 

S{k, m) ^ l — z"-''-! (E.ll) 
Jc 27ri 

where both k and m are assumed to be positive integers and the integration in the complex plane goes over a closed 
contour C around zero. Using this representation we obtain 



M 00 



27ri ^-^ V z 

Q — 1 ric — 1 



where it is assumed that the radius Rc of the the contour C is big enough, namely Rc > maxc, \aa\, so that all the 
complex numbers ai,a2,...,aM are contained inside C. In this case all the summations in Ea. (|E.12[) are convergent, 
and we find 

SiN,M) ^ i ^ z^-i . . ... . (E.13) 

Jc 2711 (z-ai) (z-a2) [z - aM) 

Since the above integral is equal to the sum of M pole contributions at z = (a = 1, 2, ...M), one gets Ea. (jE.10p . 
Now substituting Ea. (jE.10p (with = exp[A(j/Q + T-Va)]) into Ea. (IE.9|) we obtain 



Z(iV,A) - J V,{y,p) exp(AiV2/) + 

M=2 ' Q=2 ~ ^ 

where we have used the symmetry of the expressions in Eqs. (|E.7[) . (jE.8|) with respect to permutations of {{ya + iija)}- 



Appendix F 



In this Appendix, we show that, independently of the values of the momenta Pa and the parameter s, the integro- 
differential factor in Eq. (j87p provides unity, 

(F.l) 
(F.2) 



Substituting here the definitions, Eqs (|E.6|) and (|E.7I) . the factor Q can be split into two factors 

G — Ga Gb 



where 



and 



Gb = 



M 

n 



■ M 

n 

a^2 



+00 



27r dxaf} 



-\-oo 



dlpald(j)al d 



+00 



d^piadcpia d 



X exp 



M 



27r dxia 

M 



27r dxai 



Q = 2 



a=2 



A/ 



M 



{Xal,Xla}=0 



(F.3) 



(F.4) 
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For the factor Qa, Ea. (IF.3[) . we easily find 

M 



n 



1 



In a similar way we obtain for the factor Q b 

M ^+00 



\Pa-Pp\ = 1 



(F.5) 



n/ ^exp(-i[K-p,i-..]^L)xn/ 



a=2 
M 



^"^^^ exp( --^[ba -Pil +is](j}la ) X 



TT ^ exp ( [\pa-pi\- is] Xal 



d 



X ^ exp( [\p„ -pi\ +is]xic 

Xol=0 a=2 1" 



Xio=0 



fr i , 



X [(Pa -Pl)^ +s2] = 1 



(F.6) 



(F.7) 



which proves Ea. (|F.lll . 
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